Nucleation and Growth of GaN/AIN Quantum Dots 
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We study the nucleation of GaN islands grown by plasma-assisted molecular- beam epitaxy on 
AlN(OOOl) in a Stranski-Krastanov mode. In particular, we assess the variation of their height and 
density as a function of GaN coverage. We show that the GaN growth passes four stages: initially, 
the growth is layer-by-layer; subsequently, two-dimensional precursor islands form, which transform 
into genuine three-dimensional islands. During the latter stage, island height and density increase 
with GaN coverage until the density saturates. During further GaN growth, the density remains 
constant and a bimodal height distribution appears. The variation of island height and density 
as a function of substrate temperature is discussed in the framework of an equilibrium model for 
Stranski-Krastanov growth. 
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I. INTRODUCTION 

Zero-dimensional semiconductor quantum dots (QDs) 
have attracted much interest in the last decade due 
to their multiple potential applications ranging from 
low-threshold laserai*2, via single-electron tunneling 
devices^ to possible realizations of qubits for quantum 
computation!^ A versatile method for the fabrication 
of semiconductor QDs is their self-assembled growth fol- 
lowing the Stranski-Krastanov (SK) growth modei This 
mode usually occurs during the growth of semiconductor 
epilayers under compressive strain. Example material 
systems are In^Gai^As/GaAs^SuUii Si^Gei-^/SitiiiiS 
CdSe/ZnSe^ or GaN/AlNiiiSi In this mode, atoms are 
initially deposited in form of a two-dimensional pseudo- 
morphic wetting layer. The associated strain energy in- 
creases with the thickness of the wetting layer and is 
finally clastically relieved by the formation of islands. 

The usefulness of such self- assembled nanostructures 
relies on the ability to obtain homogeneous size distribu- 
tions as well as to control their size, density and position. 
Many theoretical contributions have enhanced our under- 
standing of the size distributions of SK-grown islands, 
but some controversy remainsii£iiLi£ii2i22i2ii2£ Also, de- 
spite many experimental studies, the influence of growth 
parameters on the size and density of such islands is not 
fully understood, owing to the complexity of the physics 
of strained layer growth. 

In this work, we present results on the nucleation of 
GaN islands on A1N following an SK mode, in particular 
on the dependence on the amount of deposited GaN and 
the substrate temperature. We find that the qualitative 
behavior is similar to that found in other systems, e.g. 
for InAs/GaAs and Ge/Si. This reinforces the idea that 
there are common features of semiconductor SK growth, 



which are rather universal and independent of the spe- 
cific material system. Absolute island sizes then densi- 
ties will depend on the material system, possibly through 
material parameters like lattice misfit, elastic constants, 
or surface energies. Hence, we may further discuss the 
experimental data in the framework of an equilibrium 
model based on statistical physics that was originally de- 
veloped in work on the Ge/Si system^ 



II. EXPERIMENTAL 

The samples have been grown in a MECA2000 
molecular-beam epitaxy (MBE) chamber equipped with 
standard effusion cells for Ga and Al evaporation. The 
chamber also contains an rf plasma cell to provide ac- 
tive nitrogen for GaN and A1N growth. The pseudo- 
substrates used were about 2 /jm thick GaN(OOOl) (Ga- 
polarity) layers grown by metal-organic chemical vapor 
deposition on sapphire. The substrate temperature T$ 
was measured by a thermocouple in mechanical contact 
to the backside of the molybdenum sample holder. To 
ensure substrate temperature reproducibility, each series 
of samples described below was grown on a single molyb- 
denum substrate holder. 

Prior to all experiments, a 100 nm thick GaN layer was 
grown under Ga-rich conditions on the pseudosubstrates 
to avoid the influence of a possible surface contamina- 
tion layer. Subsequently, a 300 nm thick A1N film was 
deposited under Al-rich conditions at a substrate tem- 
perature of 730 °C. We have found by reflection high- 
energy electron diffraction (RHEED) and high-resolution 
X-Ray diffraction that this thickness is sufficient for the 
A1N layer to be virtually fully relaxed with a residual in- 
plane strain of ei < 0.1%i24 An atomic-force microscopy 
(AFM) image of such an A1N(0001) surface is shown 
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FIG. 1: Atomic- force micrograph of an A1N pseudosubstrate. 



in Fig. n The surface is characterized by about 30 run 
wide terraces and spiral hillocks, similar to GaN surfaces 
grown under equivalent conditions. 25 

The growth rate and the GaN coverage have been ex- 
perimentally determined for each sample by RHEED os- 
cillations occurring during the growth of the wetting layer 
prior to island formationiSi Typically, growth rates for 
different layers (and different samples) were reproducible 
within about 0.01 monolayers (ML)/s. The GaN coverage 
was then calculated with a precision better than 0.1 ML. 
It is worth noting that no measurable influence of the 
growth rate on GaN island properties has been found in 
the range between 0.1 and 0.6ML/s. As, on the other 
hand, the Ga/N flux ratio has been found to have a cru- 
cial influence on the growth mode^ it was fixed to 0.8 
(N-rich conditions) , which leads to a critical thickness for 
the SK transition of 2.25 MLfSL 2 ! The critical thickness 
was measured for each sample and found to be repro- 
ducible within 0.1 ML. 

To study the GaN islands as grown, the samples were 
rapidly quenched to room temperature under an N-fiux. 
Ex-situ AFM was then used to study the GaN morphol- 
ogy after exposure of the samples to air. 
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FIG. 2: Atomic-force micrograph of the GaN surface mor- 
phology for coverages of (a) 1.8 ML and (b) 2.2 ML, re- 
spectively. In (b) , d ss 20 nm indicates the typical lateral 
length scale of the 2D precursor islands. 



FIG. 3: Atomic-force micrograph of GaN surface morphology 
for coverages 9 of (a) 2.5 ML, (b) 2.8 ML, (c) 3.0 ML, and 
(d) 4.6 ML, respectively. 



III. RESULTS 

A. Dependence on GaN coverage 

The parameter most directly governing the properties 
of GaN islands grown in an SK mode on A1N is the GaN 
coverage 0, which effectively describes the time evolu- 
tion of the islands during growth. To study this evolu- 
tion, a series of samples has been grown at a substrate 
temperature of 730 °C and a growth rate of 0.15 ML/s. 
The GaN coverage was varied between 6 = 1.8 ML and 
9 = 4.6 ML. 

Figure Ufa) shows an AFM image of the morphology 
of GaN layers obtained after the deposition of 1.8 ML, 
i.e. for a coverage well below the critical thickness of 
2.25 ML. We find that the morphology is unchanged with 
respect to that of the A1N pseudosubstrate: the surface 
is characterized by about 30 nm wide terraces and spiral 
hillocks. We can thus infer that the growth of about the 
first 2 ML of GaN occurs in a layer-by-layer mode since 
RHEED oscillations are observed^ 

When the second monolayer is completed, the mor- 
phology changes, as evidenced in Fig. Efb) for = 
2.2 ML, i.e. immediately before the SK transition. Re- 
mainders of terraces and spiral hillocks are still visi- 
ble, but the surface is characterized on a short scale by 
1-2 ML high flat 2D islands with typical diameters of 
d ~ 20 nm (see Fig. Efb)). The behavior can thus be 
described by a transition from a layer-by-layer growth to 
multilayer growth at around 2.0 ML. Such an occurrence 
of 2D precursor islands prior to the genuine 2D-3D SK 
transition has also been observed in the InAs/ (Al,Ga)As 
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FIG. 4: Height distribution of GaN islands grown on A1N 
with nominal coverages as indicated. Ts = 730 °C. 
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FIG. 5: Average island height as a function of GaN coverage 
at Ts — 730 °C. Above 2.8 ML a bimodal size distribution is 
observed. 



the density of mode 1 islands is identical with the to- 
tal density and increases strongly with coverage. After 
the deposition of 2.8 ML, the density of mode 2 islands 
increases strongly, similarly to the behavior of mode 1 
islands after the 2D-3D transition, whilst the density of 
mode 1 islands decreases. As the total island density 
remains approximately constant, mode 1 islands trans- 
form into mode 2 islands, probably without additional 
nucleation of new (mode 1) islands. 



system, 28 ! 29 ! 30 ! 31 

When the growth is continued, the genuine 2D-3D 
transition occurs and 3D GaN islands are formed as 
shown in Fig. [21 We find that their density increases 
strongly with GaN coverage and saturates around 3.0 ML 
at a value of 1.8 x 10 11 cm~ 2 [see Figs.EJa) and (b)]. Fur- 
ther GaN deposition does not lead to an increase in island 
density but instead islands grow in size. However, the is- 
lands do not grow continuously in size but a bimodal size 
distribution is observed [see Figs. E3[c) and (d)]. 

The variation of the island size is summarized in Figs.Q] 
and|S] We see that, initially, the islands' height increases 
from 0.7 nm to 1.6 nm between 2.4 and 2.8 ML of GaN 
coverage. For higher GaN coverage, a bimodal size dis- 
tribution is observed with the average height of the first 
mode islands remaining constant at 1.6 nm independent 
of coverage. In contrast, the average height correspond- 
ing to the second mode increases continuously reaching 
4.2 nm at 4.7 ML coverage with no sign of saturation in 
the examined coverage range. Differences in the shapes 
of mode 1 and mode 2 islands will be discussed below. 

The variation of the island density is depicted in Fig. [5] 
We observe that the total density increases strongly after 
the 2D-3D transition but saturates after 2.8 ML, i.e. at 
the coverage where the island size distribution becomes 
bimodal. The partial densities of the two modes are 
shown in the inset in Fig. EI For coverages below 2.8 ML, 



B. Dependence on substrate temperature 

The influence of the substrate temperature is stud- 
ied in a series of samples with 9 = 3.0 ML deposited 
at substrate temperatures between Ts = 690 °C and 
Ts = 760 °C. At lower temperatures, no SK growth is ob- 
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FIG. 6: Total GaN island density as a function of GaN cov- 
erage at Ts — 730 °C. The inset shows the partial density of 
the islands in the two modes. We observe that mode 1 islands 
transform into mode 2 islands after the deposition of about 
2.8 - 3.0ML. 
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FIG. 7: Atomic-force micrograph of GaN surface morphology 
as a function of substrate temperature: (a) Ts = 700 °C; (a) 
T s = 760 °C 



served and GaN grows in a pseudo-2D modeji^2£i2! prob- 
ably due to excessively low adatom mobility^ Higher 
substrate temperatures are prohibited by our experimen- 
tal apparatus due to indium bonding of the substrates 
to the substrate holder. All samples in this series were 
grown at a growth rate of 0.25 ML/s. 

An AFM image of the morphology of two GaN lay- 
ers grown at substrate temperatures of Ts = 700 °C and 
Ts = 760 °C, respectively, are shown in Fig. We ob- 
serve that the total island density decreases rapidly with 
substrate temperature, in keeping with previous results^ 
and also with results obtained for the InAs/GaAsS and 
Ge/Si^a systems. The variation of the GaN island den- 
sity as a function of substrate temperature for a coverage 
of 9 = 3.0 ML is depicted in Fig. |SJ We find an approx- 
imately exponential decrease of the island density with 
increasing substrate temperature from 4.1 x 10 11 cm~ 2 at 
T s = 690 °C to 8.0 x 10 10 c m - 2 at T s = 760 °C. 

As we have seen above, the island density tends to 
saturate at sufficiently high GaN coverage. This satura- 
tion density is shown as a function of substrate temper- 
ature in the inset in Fig. |S| For substrate temperatures 
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FIG. 8: GaN island density p for a coverage of <d — 3.0 ML 
as a function of substrate temperature Ts. The inset shows 
the island saturation density as a function of substrate tem- 
perature. 
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FIG. 9: Height distribution of GaN islands grown on A1N at 
substrate temperatures as indicated. O = 3.0 ML. 
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FIG. 10: Average island height as a function of substrate 
temperature for a GaN coverage of = 3.0 ML. Above 720 °C 
a bimodal size distribution is observed. 



of Ts < 730 °C, the saturation density is similar to the 
density for = 3.0 ML, in agreement with the results in 
the last section. However, for higher substrate temper- 
ature, the saturation density is significantly higher, in- 
dicating that saturation occurs for GaN coverages larger 
than 3.0 ML. 

Figure shows the island height distribution for = 



5 



3.0 



£ 2.0 
CD 



JO 
CO 



1.0 



0.0 







I I 




- 


• = 2.5 ML 








■ = 3.0 ML 




EED 






















i I i i i i I 


(a) : 



5 10 15 

Island diameter (nm) 



20 



D) 
CD 

sz 
"O 

c 
jo 

CO 




10 15 20 25 30 

Island diameter (nm) 



35 



FIG. 11: (a) Island height as a function of diameter for 
mode 1 islands obtained for two GaN coverages, as indicated. 
Ts — 730 °C. The solid line corresponds to an aspect ratio 
of t = 0.12. (b) Island height as a function of diameter for 
mode 1 and mode 2 islands. T s = 730 °C, 9 = 4.0 ML. The 
solid lines correspond to aspect ratios of r = 0.13 for mode 1 
and r = 0.19 for mode 2, respectively. 



3.0 ML and substrate temperatures as indicated. The 
behavior is summarized in Fig. I1UI We find a single ap- 
proximately Gaussian distribution (in the limit of the 
statistical precision) at a height of 1.6 nm for substrate 
temperatures Ts < 720 °C. These islands correspond to 
the mode 1 islands observed at Ts = 730 °C and dis- 
cussed above. At low temperature, the distribution is 
thus still monomodal after 3.0 ML. For T s = 730 °C, a 
shoulder appears at the high island side of the distribu- 
tion and transforms into a clearly separated second mode 
at higher temperatures. Hence, we find that bimodal 
distributions occur earlier (for lower GaN coverage) at 
higher substrate temperature. Another remarkable find- 
ing is that the height of the mode 1 islands remains con- 
stant as a function of temperature, whereas the height of 
mode 2 islands increases with temperature, again without 
signs of saturation in the examined temperature range. 



C. Discussion 

The above results demonstrate the occurrence of bi- 
modal island size distributions at high GaN coverage 
and/or substrate temperature; at low GaN coverage 
and/or substrate temperature, the size distribution is 
monomodal. Remarkably, when bimodal distributions 
occur, the size of mode 1 islands (smaller islands) ap- 
pears independent of growth parameters such as GaN 
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FIG. 12: Total volume contained in 3D GaN islands as a 
function of coverage. The solid line indicates the expected 
behavior for a constant 2 ML thick wetting layer, whereas the 
dashed line indicates the expected behavior for no wetting 
layer at all. 



coverage and substrate temperature. 

As discussed above, the variation of partial densities of 
both modes implies that first mode 1 islands are nucle- 
ated, their density saturates and the growth proceeds fur- 
ther by transformation of mode 1 islands into mode 2 is- 
lands. Thus the question arises, whether this transforma- 
tion involves a shape change of the islands similar to the 
behavior observed in the Ge/Si system, where pyramidal 
islands are observed to transform into "dome" -shaped 
islands during growth^S Some information about the 
island shape can be gathered by RHEED. The RHEED 
images corresponding to pure mode 1 island morphologies 
and morphologies with a mixture of mode 1 and mode 
2 islands (but with a majority of mode 2 islands) are 
both characterized by the same {1013} facets, showing 
a six-fold rotation symmetry. It thus appears that both 
types of islands are (truncated) hexagonal pyramids with 
{1013} sidewalks, in agreement with previous resultsii^i^ 
However, it is clear that the shape of such islands cannot 
be obtained from the RHEED pattern alone and further 
work is necessary to determine the precise shape of both 
types of islands. 

A more detailed analysis can be done by extracting 
the aspect ratio r of the islands from the AFM data. 
Figure lLTT a) depicts the aspect ratio of mode 1 islands for 
two different GaN coverages (at which monomodal height 
distributions are observed) at a substrate temperature of 
730 °C. In spite of the absolute increase of the islands' 
size, their aspect ratio remains constant and we find r 
0.1 '1. As shown in Fig. Illf bL the aspect ratio of mode 1 
islands remains constant (r = 0.13) after the transition 
to a bimodal distribution but the mode 2 islands have a 
significantly larger aspect ratio of r = 0.19. Both aspect 
ratios correspond to truncated pyramids^ but mode 1 
islands are flatter. 

From the aspect ratio data, we can further calculate 
the amount of GaN contained in the 3D islands as a func- 
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tion of deposited GaN. The result is shown in Fig.^J We 
find that the data are consistent with a 2 ML thick wet- 
ting layer, independent of the amount of deposited GaN. 
It thus appears that in this system, the wetting layer is 
stable and does not contribute to island growth. The 
deviation for GaN coverages just after the critical thick- 
ness of 2.25 ML can be explained by the presence of 2D 
precursor islands, which still contain significant material 
but are not taken into account as 3D islands. 



IV. MODELING 

In order to understand more precisely why the ob- 
served island distributions form, we analyze the nitride 
nanostructure growth in terms of a thermodynamic equi- 
librium model that has so far only been applied exten- 
sively to the Ge/Si system^ To date, little effort has 
been made to fit data on size distributions of nitride 
nanostructures to any model. Nakajima et alSL have 
undertaken a purely theoretical study and explicitly de- 
rived thickness-composition phase diagrams for the ex- 
pected growth mode for GalnN/GaN and GalnN/AIN, 
using known materials parameters. However they did 
not take into account the shapes and sizes of experimen- 
tally observed nitride nanostructures, nor did they seek 
to predict size distributions. Here we take a rather differ- 
ent approach. We apply the thermodynamic model to the 
AFM data on nitride nanostructure distributions grown 
at one particular temperature in order to assess whether 
the bimodal distributions of nanostructures are consis- 
tent with thermodynamic equilibrium, as opposed to be- 
ing configurations along an unstable ripening trajectory. 
The experimental distributions are described reasonably 
well by the model, as explained below. We then apply the 
model to understand the variations in the bimodal distri- 
butions with growth temperature at fixed coverage, and 
again find reasonable agreement over a range of growth 
temperatures. 

The thermodynamic equilibrium model makes explicit 
predictions for the form of the bimodal epitaxial nanos- 
tructure distributions. The model is governed by a set 
of parameters that describes the size dependence of the 
internal energy of an epitaxial nanostructure. Following 
Shchukin et al.r&we describe the internal energy ei of the 
ith individual island of type X and volume z/, as follows: 

e i = A x u i + B x ^ /3 + C x ^ /3 + D x . (1) 

Using this expression for the internal energy, the island 
size distribution f{vi) is given by 

where fj, is the chemical potential and ks is Boltzmann's 
constant. 



Equation JIJ describes the chemical and elastic con- 
tributions to the island energy. The first term arises 
from bulk strain, the second from surface and interfa- 
cial energies and island-island interactions, and the third 
from surface stress and edge energies (suppressing the log 
dependence of the edge elastic relaxation energy). The 
equation also includes a size-independent term that is 
relevant only insofar as it differs from one nanostructure 
type to another. For the internal energy JQ), a mini- 
mum may exist if B and C have opposite signs. There 
has been some debate over which of the two constants is 
likely to be positive and which negative. In the initial 
development of a "shape- map" for the bimodal Ge/Si 
system by Williams et al., ls B was assumed to be nega- 
tive and C positive. In more recent work on Ge/Si, Rudd 
et alm± have suggested that B should be positive, since a 
Ge(OOl) surface does not spontaneously roughen, though 
they also acknowledge that, since 3D growth is preceded 
by the appearance of increasing numbers of defects in the 
wetting layer with their own associated energy, the be- 
haviour of the Ge(OOl) surface of a bulk crystal may not 
be relevant and B could be negative. Here, we find that 
the option most consistent with the nanostructure size 
distributions observed by AFM is B > and C < 0. 

Unlike in the most recent developments of the equilib- 
rium models we have not, as yet, explicitly included the 
effects of the elastic interaction between pairs of islands. 
We also have not undertaken a self-consistent calculation 
of the chemical potential. Instead, we have taken a sim- 
ilar approach to Williams et al.^ Using the expressions 
for the distributions J2J) and internal energies , we have 
fitted the model to the available data at several coverages 
and temperatures, in order to determine approximately 
how the constants A', B, C, and D vary as a function of 
growth parameters. For convenience, we have introduced 
A 1 = A — [i. We then used our empirical functions to fit 
the model across multiple data sets, to determine how 
well it compares with the overall behaviour of the sys- 
tem. By comparison with the more sophisticated work 
of Rudd et alm& on Ge/Si, we expect B to be dependent 
on the coverage and A' to depend on the chemical poten- 
tial, a non-trivial function of coverage and temperature. 



A. Evolution of size distribution with coverage 

In investigating the variation of the parameters A' , B, 
C, and D as a function of coverage, we found that C and 
D could be treated as being independent of coverage, 
but that B increased monotonically with increasing cov- 
erage. This coverage dependence is consistent with elas- 
tic island-island interactions, which are expected to scale 

2/3 

with v i , and which will increase with increasing cover- 
age. B (0) could be approximated to a straight line with 
non-zero slope for both the smaller (S 1 ) and the larger 
(L) island types. Thus (for X 6 {S,L}): 
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B x (Q)~B 0X + b x e. 



(3) 



Similarly, Williams et al*& found that, for the Ge/Si 
case, C was independent of coverage and B varied lin- 
early with coverage. However, they were also able to ap- 
proximate the variation of A! with coverage as a straight 
line with non-zero slope. This does not appear to be valid 
in this case, and we have not found a simple functional 
relationship between A' and coverage. Hence in fitting 
the model across multiple data sets we have calculated a 
separate value of A' for each data set, but have assumed 
a constant relationship between the value of A' for the 
larger and smaller islands. Essentially: 

A' L (G)=g(e) and A' S (Q) = g(G) + A , (4) 

where Aq = As — A L , representing the intrinsic differ- 
ences in bulk elastic energy between the two island types. 
The function <?(0), which has not been found explicitly, 
represents the coverage dependence of the chemical po- 
tential. 

Hence, in fitting across four coverages at which large 
and small islands coexist, we used the parameters Bos, 
Bql, bs, b L , C s , C L , D s , D l , A , and four separately 
fitted values of 3(6). This clearly gives us a large num- 
ber of parameters — however, the model must fit a large 
number of observable features of the data. We are able 
to provide a reasonable fit for each shape for the rela- 
tive height, position, width and skewness of the volume 
distributions (essentially 15 experimental observables for 
each island shape). Example fits are shown in Fig. ^] 
The fitting has been done using a \ 2 approach and the 
goodness-of-fit may be judged by calculating an overall 
Q- value across the data sets, Qtotai- For the data sets un- 
der consideration, Qtotai = 0.22. Since the available data 
sets are much smaller for this system than for the Ge/Si 
system (in which counting large numbers of islands is fa- 
cilitated by the larger sizes of the islands, and the clear 
shape difference between island types) this Q-value is sat- 
isfactory. For comparison we have also attempted to fit 
the data with the original model of Williams et aiM* and 
with some Gaussian functions. These fits were markedly 
worse and it was not possible to find parameters that fit 
the entire set of data with a reasonable Q-value. From 
these findings we may infer that the data are consistent 
with the equilibrium model; however, the model can in 
no way be said to have been proven. The model contains 
a large number of parameters for which little a priori 
information is available in the nitride systems, and the 
measured island distributions would benefit from more 
data in order to reduce the statistical uncertainty. Hence, 
the discussion that follows is somewhat tentative. 

Using the calculated parameters we can plot the inter- 
nal free energy versus volume curves for both the smaller 
and the larger islands. The results are shown in Fig. 1141 
for three different coverages. There is a significant evo- 
lution of the curves as the coverage increases. Similar 
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FIG. 13: Comparison of experimentally observed island size 
distributions at 730 °C with fitted island size distributions 
for both the smaller (circles) and larger (crosses) islands, for 
coverages of 3.0 ML and 3.7 ML, as indicated. The fitting 
has been done, as described in the text, with a single set of 
parameters to describe a range of coverages. Values of the C- 
and D-parameters do not vary with coverage, values of the 
B-parameters vary linearly with coverage and values of the 
A-parameter are individually fitted such that As — Al = Aq 
for all coverages. 



results were observed by Williams et ali& for the evolu- 
tion of internal energy with coverage in the Ge/Si system, 
and were attributed to the interactions of the nanocrys- 
tals. The value of B(Q) affects the width of the curve 
and the horizontal position of the energy minimum. In 
the Ge/Si case, Williams et al. found that for the larger 
islands (domes) the value of J5(9) increases (becomes 
less negative) with increasing coverage, whilst the value 
of B(Q) for the smaller islands (pyramids) decreases, 
which suggested that the repulsion between domes was 
stronger than the repulsion between pyramids. In con- 
trast, the work of Rudd et a/,— (also on Ge/Si) suggests 
that B(Q) increases with increasing coverage at a similar 
rate for both pyramids and domes. Similarly, in the cur- 
rent case, we find that the i?-parameter increases with 
increasing coverage at almost the same rate for both is- 
land types, indicating no significant difference between 
the elastic island-island interactions for each island type. 

Changes in A' (0) tend to shift the internal- free en- 
ergy versus volume curves vertically. At low coverages, 
the minimum internal energy is lower for the smaller is- 
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FIG. 14: Plots of the variation in internal energy of an is- 
land with island volume for the smaller (solid line) and larger 
(dashed line) islands at 730 °C for three coverages, as indi- 
cated. 



lands, but at high coverages, it is lower for the larger 
islands. Hence the formation of large islands rather than 
small islands is preferred at high coverages. However, 
the elastic interactions between islands may mean that 
not all small islands grow to become large islands. Some 
may be destabilised by neighbouring islands and dissolve. 
One surprising feature of the fits is that A' is found to 
be positive for the larger islands, for all coverages. Since 
the bulk elastic energy is expected to decrease when an 
island is formed, this suggests that the chemical poten- 
tial is negative and relatively large (since A' = A — //). 
This finding is not yet well understood. 



B. Evolution of size distribution with temperature 

In order to understand the variation of the size distri- 
bution with temperature, we must consider the temper- 
ature dependence of the free energy of each island type. 
In principle, A', B, and C may all be temperature de- 
pendent. The size-independent D-parameter is constant 
and is found to be the same for both small and large is- 



lands. Williams et idM- performed some limited studies 
on the temperature dependence of the fitting parame- 
ters in Ge/Si and found that only A' exhibited signifi- 
cant temperature dependence. Similarly, in this case, C 
is found to be essentially independent of temperature. 
However, for both island types, we have found that B 
exhibits a significant temperature dependence, such that 
for X 6 {S, L}: 



B X (T S ) = B 1X + PxT s 



(5) 



The B-parameter increases at a similar rate for the 
larger and the smaller islands, implying that the surface 
energy contribution from each island type varies in a sim- 
ilar way with temperature. This is congruent with the 
suggestion that the {1013} facets seen in RHEED dom- 
inate both island types. These results are not, however, 
incompatible with the existence of other facets on the 
larger island types. 

The A'-parameter, which was found to exhibit a com- 
plex dependence on coverage, also exhibits a complex de- 
pendence on temperature. Hence in fitting the model 
across a range of temperatures, we have calculated val- 
ues of A 1 individually for each data set, but have assumed 
a constant relationship between the value of A 1 for the 
large islands and the value of A' for the small islands, 
such that 



A' L (T S ) = h(T s ) and A' S (T S ) = h(T s ) + A ± . (6) 

For the larger islands, A' is again found to be posi- 
tive at lower temperatures but decreases smoothly with 
temperature and changes sign at « 755 °C. The fits are 
shown in Fig. ^3 Qtotai = 0.73. The function h(Ts), 
which has not been found explicitly, includes the tem- 
perature dependence of the chemical potential, and any 
variation in the elastic relaxation with temperature. 

This model is only valid if a quasi-equilibrium may be 
assumed. This assumption is particularly problematic 
with regards to variation in temperature, since in par- 
ticular at lower temperatures, the system may not reach 
equilibrium, whilst at higher temperatures significant in- 
tcrdiffusion between the islands and the substrate may 
occur. This alloying will effectively reduce the mismatch 
between the islands and the substrate - strongly affect- 
ing the island size distributions, since the mean equilib- 
rium nanostructure volume varies inversely as the sixth 
power of the mismatch. 38 Rudd et al.'s more sophisti- 
cated workjS in creating the nanostructure diagram for 
the Ge/Si system, could be augmented to account for 
uniformly distributed alloying, and for kinetic effects. 
In the past, whilst the temperature dependence of the 
equilibrium island size distribution has been considered 
theoreticallyji&S very little high quality data has been 
available for comparison. The experimental studies by 
Williams et ali& only involved temperatures of 550 °C, 
600 °C, and 650 °C, and there was considerable uncer- 
tainty in the temperature measurement. Additionally, 
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FIG. 15: Comparison of experimentally observed island size 
distributions with fitted island size distributions for both the 
smaller (circles) and larger (crosses) islands, for a coverage 
of 3.0 ML and temperatures 730 °C and 750 °C, as indicated. 
The fitting has been done, as described in the text, with a 
single set of parameters to describe a range of temperatures. 
Values of the C- and D-parameters do not vary with temper- 
ature, values of the _B-parameter vary linearly with temper- 
ature and values of the yl-parameter are individually fitted 
such that As — Al = A\ for all coverages. 



the sample grown at 550 °C is believed to have been 
significantly influenced by kinetic effects, and the sam- 
ple grown at 650 °C experienced significant intermixing. 
Hence the data considered here, relating to a smaller 
range of more accurately measured temperatures, are 
very valuable as a first quantitative test of the tempera- 
ture dependence of the equilibrium model. In addition, 
the GaN/AIN system presents an important advantage 
over the Ge/Si system since it has been show that alloy- 
ing is very limited even at high growth temperatures^ 
In considering the results of the fitting, we notice that 
fitting across the temperature data sets gives a higher 
Qtotai than fitting across the coverage data, with an es- 
sentially similar fitting methodology. One might thus 
speculate that, since the samples used for the coverage 
data set were grown at the relatively low temperature 
of 730 °C, full equilibrium has not been achieved and ki- 
netic effects are having some influence. A worthwhile 
approach, in investigating this further, will be to exam- 
ine a series of samples that have been annealed at growth 



temperature to allow further evolution towards equilib- 



V. CONCLUSION 

Using ex situ AFM, we have studied the evolution 
of thin GaN layers grown by plasma-assisted MBE on 
A1N(0001) layers at substrate temperatures between 700 
and 750 °C. Initially, 2 ML of GaN grow in a 2D mode, 
followed by the occurrence of 2D islands. These islands 
act as precursors for 3D islands, which appear after an 
SK transition around 2.3 ML. During further growth, in 
particular at higher temperatures, a bimodal island size 
distribution is observed. Remarkably, the size of mode 1 
islands is found to be independent of coverage and tem- 
perature, whereas the size of mode 2 islands increases 
with coverage and temperature. The analysis of the par- 
tial island densities reveals that, whilst the total island 
density remains constant, mode 1 islands transform dur- 
ing growth into mode 2 islands. The aspect ratio of the is- 
lands is measured for both types of islands and it is found 
that they are characterized by distinctively different as- 
pect ratios, whereas no additional facets are observed in 
the RHEED pattern due to mode 2 islands. 

These findings are examined in the framework of an 
equilibrium model for SK growth. We find satisfactory 
agreement with the experimental data, suggesting that, 
as more data become available, extending the equilib- 
rium approach may be helpful in understanding and tai- 
loring nitride nanostructure distributions. We have ex- 
amined the variation of each of the fitting parameters 
with growth conditions, and have considered how this 
may relate to the physics of this system. The calculated 
parameters appear to be compatible with the available 
data — for example the _B-parameter, relating to the sur- 
face energy of the islands shows a similar variation with 
growth parameters for each island type, which is unsur- 
prising if both island types are dominated by the same 
facet. However, the fits are in no way perfect, and this 
may be due to the influence of kinetic effects on the island 
distributions. Further experiments on the annealing of 
dot arrays at growth temperature might clarify to what 
extent the growth proceeds near equilibrium or kinetic 
effects are predominant. 
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